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Abstract. We present results of a method for an automatic search for H l shells in 3D data cubes and apply it to 
the Leiden-Dwingeloo H l survey of the northern Milky Way. In the 2nd Galactic quadrant, where identifications of 
structures are not substantially influenced by overlapping, we find nearly 300 structures. The Galactic distribution 
of shells has an exponential profile in the radial direction with a scale length of Ugsh = 3 kpc. In the 2 direction, 
one half of the shells are found at distances smaller than 500 pc. We also calculate the energies necessary to create 
the shells: there are several structures with energies greater than 10i?sN but only one with an energy exceeding 
lOOi^sN. Their size distribution, corrected for distance effects, is approximated by a power-law with an index 
a = 2.1. Our identifications provide a lower limit to the filling factor of shells in the outer Milky Way: / 2 D = 0.4 
and /sD = 0.05. 
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1. Introduction 

The interstellar medium (ISM) in galaxies is far from ho¬ 
mogeneous: it is turbulent and composed of several coex¬ 
isting phases. This creates a complex picture, so far not 
fully understood. Important constituents of this picture 
are structures known as shells, supershells and holes (usu¬ 
ally called this in association with observations in Hi ) 
or bubbles and superbubbles (in connection with Hq, or 
other wavelengths). These structures are known to ex¬ 
ist in the Milky Way Galaxy and many external galax¬ 
ies (for a summary of observations of shells see Walter 
& Brinks, 1999). Most probably they are created by en¬ 
ergy release from massive stars (winds and supernova ex¬ 
plosions), however, alternative explanation, gamma-ray 
bursts (GRBs; Efremov et ah, 1998; Loeb & Perma, 1998) 
or high velocity cloud (HVC) infalls (Tenorio-Tagle & 
Bodenheimer, 1988), have been invoked in some cases. 
The majority of the observed shells are due to star for¬ 
mation. Structures connected with the infall of HVC are 
exceptional (Ehlerova & Palous, 1996). The origin of shells 
in connection with GRBs has not been proven (Efremov 
et ah, 1999). A connection between massive stars and 
HI shells can be tested either by observation of individual 
objects or by a comparison of statistical properties such 
as galactic or size distributions. In this paper we choose 
the second method. 
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There are many papers on shells in the Milky Way 
but the majority of them are concerned with individual 
objects. The most notable exceptions are two papers by 
Heiles (1979, 1984) which contain lists of shells in the 
northern Galaxy, the paper by Hu (1981) on northern high 
latitude shells, and more recently the paper by McClure- 
Griffiths et al. (2002) on shells in the southern Galaxy. 
Unfortunately, none of these lists are complete, but are the 
best data available for the Milky Way. In this paper we try 
to create a more complete list of shells in our Galaxy. We 
use the Leiden-Dwingeloo Hi survey and create an auto¬ 
matic algorithm to find shells. Then we compare identified 
structures with shells in external galaxies and discuss their 
Galactic and size distributions. 

The structure of the paper is as follows: we start with 
a short description of the dataset used for identification 
(Sect. 2), continue with a description of the identification 
method (Sect. 3) and other searching methods (Sect. 4). 
Then we analyse results for the 2nd galactic quadrant of 
the Milky Way (Sects. 5 and 6) and summarise our find¬ 
ings (Sect. 7). In the Appendix we try to cross-identify 
the shells found with those published previously. 

2. Data 

We use the Leiden-Dwingeloo Hi survey (Hartmann & 
Burton 1997; LDS). The survey was made in 1989- 
1993 with a 25m radio telescope in Dwingeloo (the 
Netherlands). It covers 79% of the sky with a spatial res- 
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Fig. 1. Types of Hi shells. The first row shows examples 
of structures, that can be identified, the second row shows 
structures connected to the surroundings or with incom¬ 
plete walls, which are undetectable by the code. 

olution of 0.5°, a velocity resolution of 1.03 kms“^ and an 
rms noise of 0.07 K. 

We chose this survey because of high sky coverage and 
uniformity of data. LDS has not been searched for shells 
up to now. Contrary to many H i surveys used to find 
shells, like the Canadian Galactic Plane Survey (CGPS), 
LDS covers the whole sky, not only a strip around the 
Galactic equator. Therefore we are able to study shells 
with large angular sizes {Ab > 10°) and also the distribu¬ 
tion of shells at high latitudes. 

The automatic method of shell detection that we de¬ 
scribe in this paper is quite general and can be applied 
to any data cube. It is obvious that by using a data 
cube with a higher angular resolution (such as CGPS) 
we would discover more smaller structures than by us¬ 
ing LDS. However, as it will be shown later (e.g. Fig. II, 
angular dimensions of many discovered shells exceed the 
region covered by CGPS. Consequently, they would re¬ 
main undiscovered if we did not use an all-sky survey. 
Moreover, a structure with a dimension of 500 pc (which 
is quite typical) has an angular dimension larger than 5° 
at a heliocentric distance less than 5.7 kpc. Therefore, as 
we are primarily interested in shells with dimensions of 
several 100 pc, the LDS (or any other high-coverage, low- 
resolution survey) is appropriate. If we were interested in 
smaller shells, like wind-blown bubbles, a high-resolution 
survey (such as CGPS) would be better. 

Our position inside the Milky Way makes the identi¬ 
fication more difficult compared to external galaxies be¬ 
cause of substantial overlapping of material along the line- 
of-sight, because of a linear resolution that changes with 
heliocentric distance and because of problems with dis¬ 
tance determination. 

3. Searching algorithm 

A prototypical ideal H i shell has several signatures which 
distinguish it from the background: it is a region of de¬ 
creased HI emission surrounded by a dense thin wall, it is 


expanding and its shape is spherical. The reality is more 
complex: the ISM is stratified and turbulent; Hi shells 
evolving in this medium are elongated and irregular, their 
expansion slows down during the evolution; and there are 
usually intervening material, gradients induced by galac¬ 
tic rotation, and the finite resolution of the telescope, all 
of which render the identification of shells more difficult. 
Taking the difficulties mentioned above into consideration 
we conclude that the most persistent quality of an H i shell 
is the existence of a region with decreased Hi emission, 
e.g. an HI hole. Our automatic algorithm is based on this 
conclusion and it searches for holes. 

The identification procedure has three steps: I) a 
search in separate velocity channels, 2) a merging of struc¬ 
tures in adjacent velocity channels to form an H i shell can¬ 
didate, and 3) a check using the central pixel spectrum. 

3.1. The search for holes in separate velocity channels 

We search for regions of local low brightness temperature 
completely encircled by higher temperature surroundings. 
Not all pixels inside the encircled region must belong to 
the hole, but all participating pixels must form a contin¬ 
uous structure. For each local minimum we search for the 
maximum encircling temperature, i.e. the one defining the 
maximum region that fulfills the described condition. Fig. 
n shows different types of prototypical structures; in the 
upper row there are structures detectable with our algo¬ 
rithm, in the lower row there are structures that are not 
detectable. The first case is a local depression surrounded 
by a wall, which has a higher density than the medium 
around it (Fig. ^ upper left). In some cases the wall is 
not observed and the structure is just a local depression 
in the H i gas. The second case (Fig. ^ upper center) is a 
local depression in the medium with a steep density gra¬ 
dient. This case is known as a galactic Hi worm and the 
structure is frequently perpendicular to the galactic disc 
because of the large-scale z-stratification of the H i . The 
last type is an arc blown out into a low density environ¬ 
ment (Fig. ^ upper right). 

All three described cases are detectable with our search 
algorithm. The difficulty begins when the local depression 
is directly connected to other low density regions (Fig. ^ 
lower left). The search algorithm is unable to close the 
structure. Thus, it either joins it to a larger neighbor¬ 
ing hole or does not identify it at all. Worms, which are 
open at one or both ends connecting regions where the 
Hi emission is absent (Fig.Q] lower center) are undetected 
as well. The identification also fails for fragmented arcs 
in a medium with low Hi emission (Fig. ^ lower right). 
Examples from the LDS of the cases described are given 
in Fig. 121 

The only constraints put on the hole are its dimensions: 
the hole must be greater than a minimum dimension and 
smaller than a maximum dimension. We use the range of 
dimensions (1.5°, 45°). There are no constraints on the 
shape of H i holes. 
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Fig. 2. Hi shells in the LDS corresponding to different types as shown in Fig.^ 


In this paper we search in lb channel maps. PV dia¬ 
grams can be used as well, however, many shells are open 
in velocity diagrams (i.e. only one hemisphere is visible) 
and this kind of structure is not well suited to our type of 
search. A full 3D search is time-consuming and also suffers 
from the mentioned disadvantages of PV diagrams. 


3.2. The formation of an Hi shell candidate 

After finding holes, i.e. continuous regions around local 
minima in individual velocity channels, substantially over¬ 
lapping holes in consecutive radial velocity channels are 
automatically joined together to form a 3D structure. If 
the structure joins four or more velocity channels (i.e. if 
its velocity extent is larger than 4 kms“^), it is called an 
HI shell candidate. In the majority of cases the shell ex¬ 
tends also to a few adjacent velocity channels, in which the 
angular size is less than 1.5°, and therefore it is not inden- 
tified by our algorithm in step 1 (but it can be identified 
in step 3, see below). If only four or five consecutive veloc¬ 
ity channels are involved, the velocity difference between 
parts of the HI shell candidate is less than the typical ve¬ 
locity dispersion in the ISM: ctism = 5-7 km s“^. In that 
case, we probably see a static structure, which will be dis¬ 
solved after the dissolution time tdiss ~ ''’sh/o’isM! where 
Cgh is the size of the shell. 


3.3. Verification using the central pixel spectrum 

As a third step in the identification scheme the spec¬ 
trum through the center of the structure is automatically 
checked. The coordinates of the center (4, be , vlsr) are 
mean values of all pixels forming the HI shell candidate. 
A ATb = Thole — Tbg spectrum through the center of an 
Hi shell candidate is analyzed for depressions and peaks. 
The background brightness temperature Tbg is defined in 
each velocity channel as an average emission from a strip 
around the structure with b = b^. This definition proved to 
be more satisfactory than other tested methods, e.g., the 
average emission from the surroundings of the structure. 

Structures that do not contain a clear depression in 
the ATb spectrum or for which the depression is not lo¬ 
cated in the velocity interval where the hole is visible in lb 
maps are excluded from further study. Fig.|31shows exam¬ 
ples of spectra. Structures in the upper row are identified 
in velocity channels where the ATb spectrum has a mini¬ 
mum, therefore all are Hi holes. The upper left gives the 
spectrum showing a minimum in those velocity channels 
where the HI shell candidate is identified (thick dashed 
line). The upper center shows a more complex spectrum 
with several smaller peaks in velocity channels correspond¬ 
ing to a temperature depression of the HI shell candidate. 
They probably come from small cloudlets engulfed by the 
HI hole. The upper right gives an example of a spectrum 
with a wall, which closes the HI hole at the position of 
the central pixel. On the other hand the three spectra in 
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Fig. 3. Spectra through the central pixel of H i shell can¬ 
didates. The upper row shows spectra where the presence 
of an Hi hole is confirmed. The lower row shows spec¬ 
tra where the signature of an Hi hole is not clear or not 
present. Solid lines are observed Tb spectra, dashed lines 
are ATb spectra, a dotted line shows the background level. 
Velocity channels in which the H i shell candidate (steps 1 
and 2) are detected are plotted by the thick line. Numbers 
are coordinates of structures (Ibv). 

the lower row do not clearly demonstrate an Hi structure. 
The spectrum in the velocity channels where the HI shell 
candidate resides is either too flat (left) or the identified 
structure is displaced from the minimum in ATb (center) 
or the structure is not connected to the minimum in ATb 
( right). 

The interval Av between peaks surrounding the de¬ 
pression in ATb or between endpoints of upward slopes 
if peaks do not exist, equals twice the expansion velocity 
of the structure Vexp = ^Av. The Av interval is usually 
larger than the velocity interval where the hole is visible in 
lb maps. This is because the algorithm searches for holes 
(i.e. minima in ATb) and moreover for holes with dimen¬ 
sions larger than a given size. Therefore, velocity channels 
that contain low-contrast or small parts of the structure 
are not included in an Hi shell candidate created in the 
first two steps of the identification. Only a minority of 
structures has an interval Av from the spectrum analysis 
smaller than the visibility interval. This is usually because 
several substructures exist inside the hole, e.g. interacting 
shells, or preexisting H i clouds. In such a case it is diffi¬ 
cult to decide which peak or upward slope corresponds to 
the rim of the structure. 

3.4. Tests on artificial data 

Tests with artificial data show that completeness in iden¬ 
tification of high-contrast shells is good regardless of 
gradients in the background. Low-contrast shells in the 
rapidly changing background can be misidentified as being 
smaller, or are completely missed. Small shells (compared 
to the thickness of the disc) are more easily identified than 
large ones. Random brightness temperature fluctuations 


create structures with a small velocity extent, which are 
easily removed when checking the reality of the feature by 
analyzing a spectrum taken through the central pixel. 

When extrapolating these results from artificial data to 
types of real shells we expect that the completeness of our 
identification for younger (i.e. not blown-out and not very 
distorted) shells is quite high, while older shells, which are 
blown-out and have fragmented walls, can be missed. We 
cannot detect shells that are formed by incomplete arcs: 
we detect only structures with a hole inside or completely 
encircled regions. 

4. Other automatic searches for H i shells 

Another approach in the search for Hi shells uses mod¬ 
els of expanding structures (usually hydrodynamic sim¬ 
ulations) and then it searches in the observed data 
cubes for patterns similar to artificial data from mod¬ 
els. This method was applied by Thilker et al. (1998) and 
Mashchenko et al. (1999) for Hi shells in the galaxy NGC 
2403, or by Mashchenko & St-Louis (2002), who use a sim¬ 
ple analytical model to search for H i counterparts of a few 
Hii regions in CGPS. Calculated shapes of Hi shells are 
usually quite smooth and regular, due to a smooth density 
distribution in the model. More complicated structures 
with irregular walls and noisy velocity patterns, which are 
results of the turbulent nature of the ISM and of preex¬ 
isting HI structures, remain undiscovered by this method. 
Our approach does not use any precomputed shapes and it 
can discover rather complicated forms. The only limitation 
is the angular and velocity resolution and the sensitivity 
limit. 

The second approach primarily uses dynamical charac¬ 
teristics of expanding shells (Daigle et al. 2003). It searches 
for a signature of an expansion in the velocity spectra, i.e. 
the existence of two peaks separated approximately by 
twice an assumed expansion velocity. Regions with these 
features are then tested if their appearance corresponds to 
an assumed shape of shells. This method was applied to a 
small field of GGPS data, where it found two expanding 
bubbles. 

In our method the spectrum through the center of an 
HI shell candidate is tested after the region is chosen as a 
depression in 16-maps, which is one difference between the 
methods. The other difference is that we do not prefer any 
expansion velocity and we do not require the existence of 
peaks in a spectrum, which is a more general condition. We 
do not put any restrictions on the shape. The spectrum 
analysis serves as a way to distinguish between random 
configurations of pixels and physical structures. 

The method as a whole is quite robust. It is not sen¬ 
sitive to the adopted model of a shell, as is the first ap¬ 
proach described above, and the spectrum through the 
structure need not contain peaks, as is the prerequisite of 
the method used by Daigle et al. (2003). Our approach re¬ 
sembles the reversed method of de Heij et al. (2002) used 
for identification of high velocity clouds in LDS. 
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The third — but from the chronological point of view 
the first — method is a traditional identihcation by eye. 
With the growing amount of data and increasing resolu¬ 
tion this method becomes more difficult and more time- 
consuming. Visual identihcation, compared to automatic 
methods, is not very effective for smaller and more regular 
or semi-regular shells. It is, however, so far unsurpassed 
when it comes to very irregular, patchy, non-expanding 
and open structures because of the ability of the human 
eye to combine disconnected features into a single shape. 

5. H I shells in the Leiden-Dwingeloo survey 

In the whole LDS we have identihed more than 600 struc¬ 
tures that passed the three steps in our identihcation 
scheme. Their positions on the sphere are shown in Fig. 
21 and their list is given in Table 1, which is available at 
the CDS. The table contains the following information: 
column 1 - the shell number, column 2 - the shell name 
in the form GSlll±bb±vvv, column 3 - the galactic lon¬ 
gitude of the center 1^ in degrees, column 4 - the galactic 
latitude of the center be in degrees, column 5 - the ra¬ 
dial velocity of the center relative to the local standard of 
rest ulsr in km s“^, column 6 - the shell dimension in 
the galactic longitude direction Al in degrees, column 7 - 
the shell dimension in the galactic latitude direction Ab 
in degrees, column 8 - the shell extent in radial velocity 
Av in km s“^, column 9 - a quality index, column 10 - 
a corresponding structure identification in other lists of 
Hi shells. Notes (1) and (2) in column 10 are explained in 
the ReadMe file of the Table. The quality index is based on 
the visual appearance of the structure and it is not part of 
the automatic search. A quality index of 1 denotes clearly 
visible and pronounced shells, while an index 4 refers to 
low-contrast and not well-pronounced shells. Indices 2 and 
3 are intermediate steps. 

There is substantial overlap of shells in the inner 
galaxy, which causes incompleteness of the identihcation. 
To avoid its inhuence on the results we decided to re¬ 
strict further analysis to 276 Hi shells discovered in the 
2nd galactic quadrant. Moreover, the 2nd quadrant is fully 
covered in LDS. 

The majority of structures have Al cosb > Ab (see Fig. 
0 ), although there is a population of ^-elongated structures 
(upward and to the left of the dotted line in Fig.[Sl). Some 
of them are worms. All large shells (> 15°) are, with one 
exception, more elongated in the I direction than in the 
b direction. This can partly be a result of the incomplete 
identihcation: because of the blowout effect the large shell 
is open to high z, and only a part of the structure is iden¬ 
tihed. This leads to an artihcally lower Ab. 

The ratio Ab/Alcos{b) is a good estimate of the shape 
of shells parallel or perpendicular to the Galactic equator. 
This ratio is not a reliable parameter for very irregular 
shells or shells elongated in directions other than the I or 
b. 

We calculated kinematic distances d using the 
Wouterloot at al. (1990) rotation curve. Some shells have 
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Fig. 5. Angular dimensions Ab against Al cosb of shells 
identihed in the 2nd Galactic quadrant. The dotted line 
corresponds to Ab — Al cosb. The upper solid line gives 
the average values of Ab/Al cosb in Al cosb bins. 
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Fig. 6. Radius of the shells Vsh versus their kinematic dis¬ 
tances. The two solid lines give the minimum and max¬ 
imum Tsh corresponding to 1.5°/2 and 45°/2 at the dis¬ 
tance d. Dashed, dotted and dash-dotted lines give the 
average Csh calculated for different models of size distri¬ 
bution of shells given by eq.(2). The histogram gives the 
average radius Cgh of the observed shells in separate dis¬ 
tance bins. 

forbidden velocities or lie close to the anticenter direction. 
This reduces the number of shells in the 2nd quadrant 
with kinematic distances to 168. The angular dimensions 
Al cosb, Ab are transformed to linear sizes. The shell ra¬ 
dius Tsh is given as 

Tsh = -^d [tan{Al cosb) + tan{Ab)]. (1) 

The dimensions of structures as a function of heliocen¬ 
tric distance are shown in Fig. 21 The solid lines indicate 
the limits imposed on Csh by restricting the range of an¬ 
gular sizes of the shells to the interval (1.5, 45) degrees. 
Let us assume that the size distribution of shells has a 
universal power-law form (Oey & Glarke 1997): 

dA^(rsh) oc rsh"“drsh, (2) 
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Fig. 4. Positions of identified Hi shells in the LDS. The cross-hatched region around I = 300° is a region not accessible 
to the Dwingeloo telescope. 


where dN is the number of shells with radius Cgh G 
{i"sh, rsh + drsh)- For shells with the size distribution ||21l we 
calculate the expected average size for a given heliocentric 
distance d as 

r”“r-“+idr 

^ min 

Integration ranges r^nin and fmax are given either by 
real physical limits (minimum and maximum possi¬ 
ble radius of the shell) or by the resolution limits 
{d tanl.5°/2,d tan45°/2) of the searching algorithm, 
whichever is appropriate (larger or smaller). The ma¬ 
jor part of the theoretical curves in Fig. 0 uses rmin = 
d X tanl.5°/2 and rmax = 1.3 kpc (if not said otherwise; 
for the choice of 1.3 kpc see Oey & Clarke, 1997). The 
thick dashed line in Fig. O refers to the size distribution 
with 0 = 2; the thick dot-dashed line to a = 3. Thin 
dotted lines have a = 2 and r„iax = 2.6 kpc (the upper 
line) or r^ax = 0.65 kpc (the lower line). According to 
a comparison with the histogram in Fig. 0 showing the 
average observed shell radii in separate distance bins, it 
seems that the power-law with a = 2 and the maximum 
radius of the shell Cgh = 1.3 kpc best fits the observations, 
but the spread in Cgh is large. 

The distribution of expansion velocities Uexp = is 
shown in Fig. [3 The largest observed expansion velocity is 
25 kms“^; the median value lies between 5 and 10 kms“^. 
There are only a few structures with Uexp < 5 kms“^, 
which is apparently because such non-expanding struc¬ 
tures cannot survive for very long in a turbulent ISM (see 
section 3.2). 



Fig. 7. Expansion velocities of identified Hi shells in the 
2nd Galactic quadrant. 

5.1. Radius vs expansion velocity 

A relation between the shell radius rgh and the expansion 
velocity ?;exp is shown in Fig. 13 We also plot lines of con¬ 
stant expansion time and luminosity over the ISM density 
{L/n) according to an analytical solution of Weaver et 
al. (1977). As expected, there are no shells younger than 
1 Myr and there are only a few shells with ages greater 
than about 50 Myr. Young shells are small and not fully 
developed, while old shells are either completely destroyed 
by external forces (like spiral arms) or too fragmented for 
our method to identify. There is a lack of small shells with 
high expansion velocities, but this is not surprising. Shells 
smaller than about 50 pc are either wind-blown bubbles 
or remnants after one supernova explosion. Such objects 
tend to have expansion velocities around or lower than 
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log r^h (^Po) 

Fig. 8. Radius of shell vs. expansion velocity for shells 
in the 2nd quadrant. Lines are analytical solutions after 
Weaver et al. (1977). Solid lines are lines of constant ex¬ 
pansion time (1, 10, 50 Myr), dashed lines are lines of 
constant L/n (0.01, 0.1, 1, 10 SN Myr“^/cm“^). Different 
symbols denote different energies of shells (according to 
Chevalier, 1974), circles are shells with E < ifsN, crosses 
are shells with E > ifsN- 

10 kms“^. Dimensions of real “supershells”, created by 
several OB stars in associations, start at radii around 100 
pc. Some supershells have expansion velocities larger than 
10 kms“^, though we do not see any clear dependence 
of radius on velocity. The argument that small shells are 
wind-blown bubbles and large shells are created by sev¬ 
eral SNs is further supported by their energies (Chevalier, 
1974). All small shells have energies lower than the typical 
energy of one supernova explosion Esn = lO^^erg. 

It is interesting to compare the properties of shells in 
the Milky Way with those in external galaxies. Walter & 
Brinks (1999; further WB99) made a comparison of four 
galaxies, M31, M33, IC 2574 and Holmberg II. They plot 
a graph of radius vs. expansion velocity (Fig. 20 in their 
paper) for these four external galaxies. They find that the 
two spiral galaxies (M31 and M33) have on average smaller 
shells and a larger range of expansion velocities compared 
to the dwarfs IC 2574 and Ho II. Neither dwarf contains 
relatively small shells (radius around 100 — 300 pc) with 
high expansion velocities (around 20 kms“^), while there 
is a significant population of such shells in M31 and M33. 
The Milky Way behaves like a spiral galaxy in this aspect; 
there are many expanding shells with a 100-300 pc size. 
Contrary to M31 and M33, in the Milky Way there are 
more shells with radii greater than 800 pc, while there is 
only one in M31 and a few in M33. This can be partly due 
to a different inclination angle of galaxies, because differ¬ 
ent viewing angles may emphasize different parts of shells. 
This effect is only important for non-spherical, therefore 
predominantly large shells. 

Differences in the rsh vs Uexp graph can be explained 
by different conditions between spiral and dwarf galaxies 
and by different star formation histories. It is possible to 


calculate the age of a shell from its size and expansion 
velocity, e.g. using the solution of Weaver et al. (1977), see 
Fig. IHl According to WB99, spiral galaxies M31 and M33 
have many young shells (younger than 10 Myr), but nearly 
no old ones (older than 30 Myr), while dwarf galaxies IC 
2574 and especially Ho H have few young shells, but many 
old ones. The mean ages of shells in these galaxies are 6 
Myr (M31), 7.2 Myr (M33), 14 Myr (IC 2574) and 37 Myr 
(Ho H); all quoted values were multiplied by a factor of 
3/5 to change them from the simple r/v used in WB99 
to the solution of Weaver et al. (1977). The differences 
are explained by different star formation histories and a 
higher rate of shell destruction by density waves in spiral 
arms. 

A somewhat different behavior was found for the 
Magellanic Clouds. Kim et al. (1999) shows an rsh vs 
'I'exp graph for the LMC. There is a population of young 
shells, but there are no shells larger than 600 pc. This is 
not caused by the missing short spacing data, as it has 
been shown by Kim et al. (2003) and Staveley-Smith et 
al. (2003): there were no new identifications of supershells 
when the Parkes and ATCA observations were combined. 
Shells larger than 600 pc are intrinsically missing in the 
LMC. The mean age of shells in the LMC is 4.9 Myr. 

This is similar to the mean age of shells in the SMC 
(Staveley-Smith et al. 1997), 5.4 Myr. The list of shells 
in the SMC also suffers from the lack of large shells, but 
this effect is not so severe: in the combined Parkes and 
ATCA data (Stanimirovic et al. 1999) only three new su¬ 
pershells with Tsh > 600 pc were found. The behavior of 
the Magellanic Clouds is clearly different from that of the 
dwarfs IC 2574 and Ho H. The presence of young shells 
in the SMC is explained by significant star formation in¬ 
duced by interaction with the Milky Way, but the lack of 
old shells - all the SMC shells are younger than 30 Myr - 
is somewhat puzzling since the last SF burst in the SMC 
happened 200 Myr ago (Stanimirovic et al. 2004). 

The mean age of shells in our list is 8.4 Myr. This value 
is influenced by different selection effects, i.e. the detec¬ 
tion of old shells is much less complete than of young ones 
since old ones are less continuous and more fragmented. 
All shells with ages greater than 30 Myr are found at galac- 
tocentric distances larger than 12 kpc — one shell with an 
age of 30 Myr is located at i? = 12 kpc, the rest of the 
old shells are at distances greater than 16 kpc. The age 
distribution of shells is similar to that for spiral galaxies 
M31 and M33. 

Our search algorithm is not very sensitive to frag¬ 
mented structures, which may explain the lack of old 
shells. However, as we do not see any old shells at dis¬ 
tances smaller than 12 kpc, while there are some at larger 
distances, we can at least say that the lifetime of shells is 
longer at greater galactocentric distances. 
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Fig. 9. An average Hi density derived for shells at differ¬ 
ent positions in the Galaxy. The grey scale is logarithmic - 
the darkest color corresponds to 0.63 cm“^, lightest color 
to 0.01 cm“^, the white color denotes positions with no 
information about the density, i.e. without an Hi shell. 
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Fig. 10. The Hi density versus size of the shells Csh- Two 
dashed auxiliary lines are power-laws with indices of -1 
and -2. 


5.2. Hi density 

The original H i density at the place of a shell is calculated 
as the density corresponding to the mass missing in the 
HI hole. If that mass were evenly added to the hole, the 
structure would disappear and the observed H i column 
density would be equal to the average value around the 
shell. The maximum value found for an individual H i shell 
is 24 cm“^ for a very small nearby shell and 2.2 cm“^ for 
shells with a radius greater than 50 pc. The typical den¬ 
sity is around 0.5 cm“^. We calculated an H i density n in 
positions of all shells and made an average for each vol¬ 
ume element of the galactic disc (Fig. EJ. As expected, 
the density decreases with increasing \z\ and R coordi¬ 
nates. These reflect the density gradients in the HI galaxy 
disk in combination with the shell size versus density ef¬ 
fect: the observations show that shells expand more in a 



Fig. 11. The distribution of energy deposited in Hi shells. 
The distribution of all shells (dotted line) and the distri¬ 
bution of shells with size Cgh >0.1 kpc (solid line). The 
distribution for the M31 galaxy (dash-dotted line) is taken 
from Walter & Brinks (1999); it is scaled to the Milky Way 
values. 

lower density medium. This is demonstrated in Fig. cni 
where we plot sizes of shells rsh versus densities n. 

5.3. Energy 


The most widely used method to estimate the energy of 
Hi shells is Chevalier’s formula (Chevalier 1974), 



where n is the density of the ambient medium, rgh the 
radius of the shell and Uexp its expansion velocity. We 
take Esn = lO^^erg. Energy estimations using Eq. 01 
of smaller and less evolved shells are quite reliable, but 
for non-spherical shells it gives an underestimation of the 
energy that created the structure. 

The distribution of the energies in H i shells according 
to the Eq. 0 is shown in Fig. IIII There are many low- 
energy structures {E < 0.1 Esn) that were probably not 
created by SN explosions. The majority of them are small 
with radii Cgh < 100 pc. Many of them are wind-blown 
shells. Even though Chevalier’s Eq. (4) has the largest de¬ 
pendence on radius, many large shells do not belong to 
the most energetic structures since the density n is sys¬ 
tematically smaller for large compared to small shells; this 
influences the energy derived with Eq. 0 in the opposite 
direction to the size Csh- 

There are only 12 H i shells with an energy E > 20 Esn 
and one with an energy exceeding 100 Esn (its energy is 
120 Esn)- This can be understood because we study the 
outer Milky Way and we do not observe many energetic 
sources there, at least not many OB associations. Bureau 
& Carignan (2002) suggest that ram pressure due to a 
galaxy moving through the inter-galactic medium of a 
group or cluster can remove material from the edges of 
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the holes and prevent it from refilling the holes’ interiors. 
Thus, ram pressure enlarges preexisting holes and lowers 
the amount of energy necessary for their creation. Ram 
pressure and high velocity cloud impacts apply particu¬ 
larly at the galactic peripheries where the gas density and 
star formation rates are low. 

We over-plot a scaled energy distribution of shells in 
M31 (a dash-dotted line in Fig. Illll taken from WB99. 
They found that the distribution is very similar in the 
four galaxies (M31, M33, IC 2574, Ho II). The Milky Way 
distribution is similar as well: it has a maximum at the 
same value (around I FIsn) and the same decline for high 
energies. However, we detect more low-energy shells than 
were found in M31, which may be partly a selection ef¬ 
fect due to the larger distance to M31, limiting the linear 
resolution. Another explanation is that the Milky Way 
creates more of them, or because our search algorithm is 
more sensitive than a visual inspection, upon which the 
M31 distribution (and all others in WB99) is based. The 
mean energy of shells in M31 is 15 Asn, similar to M33 
(16 Esn), IC 2574 (11 Asn) and Ho H (14 Asn)- The mean 
energy of shells in the outer Milky Way is 3 Asn (4.7 FIsn 
for shells with rgh > 100 pc), this lower value is related to 
the extended low-energy wing in the energy distribution 
function. 

6. Distribution of shells in the Milky Way from 

Distance Effect Removed (DER) samples 

6.1. Distance Effect Removed (DER) samples 

Due to our position inside the Milky Way our list of shells 
suffers from various distance-dependent selection effects. 
Intrinsically small structures are seen only in the vicin¬ 
ity of the Sun. When they are further out, they are too 
small to be either observed by the telescope or detected by 
our algorithm. On the other hand, large shells, when close 
to the Sun, are angularly extended and such objects are 
difficult to identify because they are too big for the field 
in which identifications take place, because of interven¬ 
ing material, and sometimes also because we usually see 
only walls of these structures and they are incomplete or 
fragmented. To avoid these problems we create ’distance 
effect removed (DER) samples’, where we include shells 
with dimensions and distances such that their angular di¬ 
mensions fall into an interval (<('min, ^max) = (1-5°; 15°). 

The lower limit dmin is the minimum dimension for 
identification, the upper limit ^max is smaller than the 
upper limit for identification (45°). This choice was made 
to ensure that we do not miss any structure in the selected 
interval. The range of heliocentric distances (or shell di¬ 
mensions) is a free parameter. When processing a DER 
sample, we work with shells whose linear dimension Agh = 
2rsh is Ash S (Amin, A„iax) and heliocentric distances d 
are d G (dmin, ^max), where Amin — dmax^a?T.(i()min) and 
Amax = dnuntan{4>ni!ix)■ With this we ensure that the 
smallest shell detected at the largest distance dmax has 
an angular size equal to ^min and at the same time that 


the largest shell at the smallest distance dmin has an angu¬ 
lar size equal to ^max- We analyze different DER samples 
selected with different choices of minimal and maximal he¬ 
liocentric distances dmin,dmax- All shells in studied DER 
samples are larger than 100 pc, which avoids contami¬ 
nation by wind-blown bubbles and low-energy structures, 
and smaller than 700 pc to avoid an incompletness in iden¬ 
tifications of large evolved shells. 

6.2. Radial distribution 

It is generally thought that H i shells are mainly produced 
by SN explosions and other energetic activities in OB as¬ 
sociations. Therefore their radial distribution should re¬ 
flect the radial distribution of young stars and OB as¬ 
sociations. Attempts to find a correlation between these 
distributions were made by Ehlerova & Palous (1996) for 
the Milky Way Galaxy and Palous & Ehlerova (1997) for 
M31 and Holmberg H. The distribution of H i shells in the 
Milky Way was constructed from the list of Heiles (1979) 

and an exponential decrease of Egsh ex exp "k®'' was 
found, with Sgsh being the surface density of Hi shells, R 
the galactocentric distance and dgsh — 4.5 kpc the radial 
scale. However, the effects of heliocentric distance were 
not taken into account properly. In this paper we will per¬ 
form a similar analysis on our new and better-processed 
list of Hi shells. 

Another attempt to determine the galactic distribu¬ 
tion of shells was made by McClure-Griffiths et al. (2002). 
Their list consists of shells from Heiles (1979) and their 
own identifications. They are mostly interested in the re¬ 
lation between shells and spiral arms. 

We fit an exponential profile to the radial distribution 
of shells identified and described in this paper: 

R 

Sgsh(-R) = EogshC (5) 

where R is the galactocentric distance and Ugsh is the ra¬ 
dial scale length. The number of shells dN at a given dis¬ 
tance R is 

dN{R) = 27rREgsh(R)di? (6) 

The observed numbers of shells in intervals of galac¬ 
tocentric distances are multiplied by a geometric factor 
proportional to the ratio of the surface of the whole galac¬ 
tocentric ring and the surface of that part of the ring that 
lies in the studied heliocentric range. These recalculated 
numbers, which correspond to dN(R) in Eq. 0, are then 
fitted to the profile 0. The resulting radial scale length 
(Tgsh for the DER sample with Amin = 200 pc and Amax 
= 800 pc is 

cTgsh = 3±i kpc. (7) 

We exclude shells with Vgh > 400 pc because DER sam¬ 
ples including these large shells would have dmin > 3 kpc. 
It would substantially reduce the number of shells in the 
interval of galactocentric distances Rqc — (9,12) kpc and 
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Fig. 12. The radial and z-distribution of H i shells. The left panel shows recomputed numbers of shells derived from a 
DER sample by multiplication with a geometric factor (see the text) at different galactocentric radii and an exponential 
fit with cTgsh = 3 kpc (solid line). The dashed line gives the numbers from the DER sample (see the text). The right 
panel is the z-distribution of H i shells. Solid lines show the normalized cumulative number of shells at coordinates z 
lower than a given value, dashed lines show the normalized energy in these shells. Binned data are observations, lines 
are exponential profiles with a scale length of 760 pc (solid) and 380 pc (dashed). Corresponding HWHMs are 530 
and 270 pc. 


increase the error of Cgsh- The fit, together with the ob¬ 
served and recalculated numbers of shells is shown in Fig. 
^l(left panel). 

The radial scale length of about 3 kpc is compara¬ 
ble to that of the stellar disc. Luminosity models give 
values between 2.5-3.0 kpc (Kent, Dame & Fazio 1991, 
(T = (3±0.5) kpc; Freudenreich 1998, a = 2.5 kpc). OB as¬ 
sociations have a radial scale length of 1.8 kpc (Bronfman 
et al. 2000). If the scale length of stars and Hi shells 
differed substantially, it would mean that an important 
part of the shells is not connected to the stellar activ¬ 
ity. This,however, is not the case. The derived radial scale 
length of HI shells roughly corresponds to the stellar scale 
length and is larger than the scale length of OB stars. The 
scale length of OB associations is based mainly on mea¬ 
surements in the inner Galaxy (for Rqc S (0.5; 2.0)Rq 
with the majority of sources inside the solar circle), while 
the scale length of HI shells is derived only for the outer 
Galaxy. Therefore it is possible that the discrepancy is 
a result of a systematic difference between the inner and 
outer Galaxy. Another possible explanation is that a cer¬ 
tain fraction of HI shells is created by a different mecha¬ 
nism than energy input from OB associations. 

6.3. The distribution in the z-direction 

Hi shells are concentrated towards the Galactic plane. 
One half of all shells lie in a 1-kpc thin layer (the HWHM 
is around 500 pc). The thickness of this layer increases 
slightly with galactocentric distance: the HWHM for shells 
with R G (9,11) kpc is 530 pc (see Fig.lT^. for those with 
R G (12,14) kpc, it is 590 pc. A more precise determina¬ 
tion of this trend is not possible because of low number 
statistics at larger galactocentric distances. 


Identified H i shells reside in a much thicker disc com¬ 
pared to the Hi disc or even to the stellar disc. It may 
be partly due to our identification procedure, which shifts 
the shell centers artificially to lower density ends, i.e. to 
higher z-coordinates. 

The average energy of shells (according to Eq. de¬ 
creases with increasing coordinate z (see Fig. 11211 . It is 
uncertain if this is a real physical effect (fewer energetic 
sources at high z coordinates) or simply a consequence of 
the fact that the volume density decreases to higher val¬ 
ues of z and as a result the estimated energy is lower. The 
HWHM of the energy release is about half of that of the 
number density, i.e. around 250 pc, also increasing with 
increasing R, 200 pc for R G (9,11) kpc and 230 pc for 
R G (12,14) kpc. 

The exponential or Gaussian fit of the z-distribution of 
HI shells is not as good as the radial one, but the exponen¬ 
tial fit is better than Gaussian. The HWHM derived from 
these fits is larger than observed. This means that there 
are more shells at high coordinates than would correspond 
to an exponential (Gaussian) profile, i.e. the distribution 
of shells has extended wings. Fig. El (right panel) shows 
the z-distribution of shells in the DER sample with Amin 
= 200 pc and Amax = 800 pc. 

6.4. Size distribution of shells 

Oey & Glarke (1997; hereafter OC97) calculate the size 
distribution of shells from theoretical considerations of 
shell evolution. They found that the size distribution 
can be quite reliably approximated by a power law @. 
They compare the size distributions of shells in the galax¬ 
ies M31, M33, Holmberg H and SMC and find that the 
power-law representation generally fits the data well with 
a e (2.1,2.7). 
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Fig. 13. The size distribution of Hi shells. The number 
of shells with radius larger than 100 pc and smaller than 
a given value of Vsh vs. the radius rsh is plotted, plus a 
theoretical value for three different power-laws. 


Fig. 14. The 2-D filling factor of Hi holes in the Tb{ 1, b, v) 
datacube for {l,b) S (90°, 180°;—10°,-1-10°) (solid line), 
{l,b) £ (90°, 180°;—5°,-1-5°) (dashed line), and {l,b) G 
(90°, 180°;—20°,-1-20°) (dotted line) as functions of the 
radial velocity. 


Fig. 1131 gives the size distribution together with three 
different power laws with a = 1.1, 2.1 and 3.0. A fit of the 
DER sample with Amin = 200 pc and Amax = 1 400 pc 
gives the average value of the index a around 2.1: 

a = (2.1 ±0.4). (8) 

The spectrum seems to be flatter for smaller shells (a < 
2.1) and steeper for larger shells (a > 2.1). There may be 
several regimes in the size distribution, one with a = 1.1, 
one with a = 2.1 and another with a — 3.0. The dif¬ 
ference in the indices may be real or it may be an ar¬ 
tifact of the identification method. Evolutionary effects 
should lead to an increase of a for large shells: some of 
the large old shells are destroyed, fragmented or other¬ 
wise distorted and therefore are more difficult to hnd. 
Merging and grouping of small shells produces a change 
in a between small and intermediate shells: instead of sev¬ 
eral small shells we see one larger shell. Our identification 
method (as discussed earlier) can underestimate the di¬ 
mensions of extended shells in the z (&) direction. Indeed, 
Fig. 0 shows that shells become more elongated in the I 
directions with increasing Al. This effect may also lead 
to higher a for large shells. An alternative explanation is 
that the completeness of our list is comparatively good for 
large shells, while it is worse for smaller shells, but this is 
unlikely as small shells are more regular and easier to find 
than large distorted structures. 

The power-law index of the luminosity distribution of 
energy sources $(£) (e.g. OB associations) is given by 
OC97 as 

$(£) oc (9) 

where C is the luminosity of the source. OC97 give the 
relation between the power laws @ and as 

( 10 ) 


which would imply from our size distribution fit 

/? = (1.6T0.3). (11) 

However, the direct connection between sizes of shells 
and powering sources disregards the size-density relation 
shown in Fig. 10. It can be one of the systematic effects, 
which distinguishes between different types of galaxies. 
The indices a and (3 are on the shallower end of the corre¬ 
sponding values given by OC97 for M31, M33, Holmberg 
H and SMC; a G (2.1, 2.7), /? G (1.6,1.9). The index /3 is 
flatter than 2, which is the slope of the H ii region lumi¬ 
nosity function (Kennicutt et al. 1989) or OB association 
luminosity function (McKee & Williams 1997). This lower 
value may be connected to the blow-out effect in the thin 
Hi disc of the Milky Way. The sizes of shells powered 
by high C OB associations are smaller because a fraction 
of the energy leaks to high-z distances. Blow-outs reduce 
the difference between the sizes of Hi shells created by 
high and low C OB associations. The shell size distribu¬ 
tion then leads to an artificially flatter luminosity function 
than that of OB associations. 

7. Filling factor of H i shells 

The filling factor / of HI shells is the total area (/ 2 d) or 
volume (/sd) occupied by HI shells compared to the total 
area or volume of the galaxy. The derived filling factors 
are lower limits to intrinsic values, since some shells are 
missing (due to distance effects, incomplete identification, 
etc). 

The 2-D filling factor as a function of radial velocity 
is shown in Fig. 1141 where we show the results for strips 
along 6 = 0°. The profiles for the wider or narrower in¬ 
tervals of b are similar. Neglecting distance effects, i.e. 
neglecting the fact that at different distances shells with 
different dimensions are not seen, the total filling factors 


a = 2/3- 1, 
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in the 2nd galactic quadrant of the Milky Way are calcu¬ 
lated. We get /sD — 0.05 and / 2 D — 0.4. We argue that 
neglecting distance effects we can approach to real, unbi¬ 
ased values. This is because filling factors are dominated 
by large shells, therefore the non-detection of small shells 
does not strongly influence the results. 

8. Summary 

We created an automatic method to search for H i shells 
in data cubes. The method is based on finding local min¬ 
ima in individual velocity channels. These 2D regions are 
then combined along the third dimension, the radial veloc¬ 
ity axis, to create 3D objects. The spectrum through the 
center of the structure is analyzed. The method has the 
advantage of being independent of models of H i shells and 
therefore being widely applicable. Small shells (compared 
to the thickness of the disc) are more easily identified than 
large ones. Random brightness temperature fluctuations 
create structures with a small velocity extent, which are 
easily discarded when looking at spectra taken through 
the centers of the shells. We cannot detect shells that are 
formed by incomplete arcs: we detect only structures with 
a hole inside or completely encircled regions. We detect 
more shells than published by Heiles (1979) and other au¬ 
thors. Some of the shells identified by us for the first time 
are as well defined as those previously catalogued. 

In the whole Leiden-Dwingeloo survey we identify 
more than 600 structures. To avoid problems connected 
with a substantial overlap of shells in the inner Galaxy, 
we analyze only structures in the 2nd Galactic Quadrant. 
There are nearly 300 shells found there. From their list 
we create samples that compensate for the distance effect, 
i.e. the fact that we see small shells only in the vicinity 
of the Sun, and that angularly large shells are difficult to 
identify in the data. 

Firstly, we compare shells in the Milky Way with shells 
in external galaxies. We find that galactic shells have sim¬ 
ilar properties to shells in M31 and M33. The majority 
of shells are ~ 10 Myr old and there are few shells older 
than ~ 30 Myr. The energy distribution is similar to that 
of M31, but we detect relatively more low-energy shells. 

Secondly, we study the Galactic distribution of 
Hi shells. The radial decrease in the surface density of 
shells in the outer parts of the galactic disk is exponen¬ 
tial, with a scale length of about 3 kpc. This is comparable 
to the radial scale length of the stellar disc, but larger than 
the radial scale length of OB associations. The similarity 
between scale lengths of stars and HI shells supports the 
idea that most shells are connected to stars. The reason 
for the discrepancy between scale lengths of H i shells and 
OB associations is not clear, it may either reflect differ¬ 
ences between the inner and outer Galaxy, or it may result 
from the fact that a fraction of shells in the outer Galaxy 
is created by a different mechanism than by an energy 
input from young and massive stars. 

HI shells are concentrated towards the Galactic plane. 
Half of all shells lie in a 1-kpc thin layer (the HWHM 


is around 500 pc). The thickness of this layer increases 
slightly with galactocentric distance. The energy in shells 
is more concentrated toward the plane; half of all the en¬ 
ergy is contained within the disc with a total thickness of 
500 pc. 

We then study the size distribution of Hi shells. We 
find that a power-law function 0 with a slope a = 
2.1 ± 0.4 describes the observed distribution, a seems to 
be larger for larger shells and smaller (flatter) for smaller 
shells. We attribute the increase of a with dimension to 
evolutionary effects and the inability of the identifica¬ 
tion method to find distorted and fragmented structures. 
The size distribution with a = 2.1 is flatter, but com¬ 
parable to the distributions found in the galaxies M31, 
M33, Holmberg H and SMG. The corresponding slope of 
the luminosity function of sources powering H i shells is 
/3 = (1.6 ± 0.3). This is also flatter than the observed dis¬ 
tribution of Hii regions or OB associations in the Milky 
Way. However, the size of shells also depends on the den¬ 
sity of the ambient medium. Gonsequently, part of the 
effect may be due to a size-density relation observed in 
our shell sample. 

The filling factors are derived for shells indentified in 
the 2nd galactic quadrant. Due to incompleteness of the 
sample, the derived values / 3 D = 0.05 and / 2 D = 0.4 are 
lower limits to the intrinsic values. 

With Tsh, Vexp, the surface density of the shell Sgh and 
with the speed of sound in the shell, it would be possible 
to estimate the gravitational stability of the shell and find 
out if it can trigger star formation. Some of the shells 
detected in the LMC, like the Sextant region near LMC IV 
(Efremov et al. 1999), probably triggered star formation. 
A discussion of shell triggering and a comparison between 
the Milky Way, LMC, Ho H and other galaxies will be 
published elsewhere. 

We conclude that our automatic identification code is 
a powerful tool in studies of Hi shells, which enabled us 
to perform the so far most extensive statistical study of 
shells in the Milky Way. 

Acknowledgements. The authors gratefully acknowledge finan¬ 
cial support by the Grant Agency of the Academy of Sciences 
of the Czech Republic under the grant No. B3003106 and sup¬ 
port by the grant project of the Academy of Sciences of the 
Czech Republic No. K1048102. They would also like to express 
their thanks to Bruce G. Elmegreen for comments on how to 
improve the English and to the anonymous referee for sugges¬ 
tions on how to improve the paper. 

References 

Bronfman, L., Casassus, S., May, J., & Nyman, L.-A. 2000, 

A&A, 358, 521 

Bureau, M., & Carignan, C. 2002, AJ, 123, 1316 
Chevalier, R. 1974, ApJ, 188, 501 

Daigle, A., Joncas, G., Parizeau, M., & Miville-Deschenes, M.- 

A. 2003, PASP, 115, 662 

de Heij, V., Braun, R., & Burton, W.B. 2002, A&A, 391, 159 
Hu, E. M. 1981, ApJ, 248, 119 


Ehlerova & Palous: H l shells in the outer Milky Way 


13 


Efremov, Yu. N., Ehlerova, S. & Palous, J. 1999, A&A 350, 
457 

Efremov, Yu. N., Elmegreen, B. G., & Hodge, P. W., 1998, 
ApJ, 501, L163 

Ehlerova, S., & Palous, J. 1996, A&A, 313, 478 
Ferriere, K. 1995, ApJ, 441, 281 
Ferriere, K. 1998, ApJ, 503, 700 
Freudenreich, H.T. 1998, ApJ, 492, 495 
Hartmann, D., & Burton, W.B. 1997, Atlas of Galactic Neutral 
Hydrogen (Cambridge University Press) 

Heiles, C. 1979, ApJ, 229, 533 

Heiles, C. 1984, ApJS, 55, 585 

Heiles, C. 1990, ApJ, 354, 483 

Heiles, C., & Habing, H. J. 1974, A&A Suppl. 14, 1 

Hu, E. 1981, ApJ, 248, 119 

Kennicutt, R.C., Edgar, B.K., & Hodge, P.W. 1989, ApJ, 337, 
761 

Kent, S.M., Dame, T.M., Fazio, G. 1991, ApJ, 378, 131 
Kim, S., Dopita, M. A., Staveley-Smith, L., & Bessell, M., S. 
1999, AJ, 118, 2797 

Loeb, A. & Perma, R., 1998, ApJ, 503, L35 
Mashchenko, S.Y., St-Louis, N. 2002, in ASP Conference 
Proceedings Vol. 260, ed. A.F.J. Moffat & N. St-Louis, 65 
Mashchenko, S.Y., Thilker, D.A., & Braun, R. 1999, A&A, 343, 
352 

McClure-Griffiths, N.M., Dickey, J.M., Gaensler, B.M., & 
Green, A.J. 2002, ApJ, 578, 176 
McKee, C.F., & Williams, J.P. 1997, ApJ, 476, 144 
Oey, M.S., Clarke, C.J. 1997, MNRAS, 289, 570 (OC97) 
Palous, J., & Ehlerova, S. 1997, in The Physics of Galactic 
Halos, ed. H. Lesch, R.-J. Dettmar, U. Mebold & R. 
Schlickeiser, 291 

Slavin, J.D., & Cox, D.P. 1993, ApJ, 417, 187 
Stanimirovic, S., Staveley-Smith, L., Dickey, J.M., Sault, R.J., 
& Snowden, S.L. 1999, MNRAS, 302, 417 
Stanimirovic, S., Staveley-Smith, L., & Jones, P. A. 2004, ApJ, 
604, 176 

Staveley-Smith, L., Sault, R.J., Hatzidimitriou, D., Kesteven, 
M.J., & McConnell, D. 1997, MNRAS, 289, 225 
Staveley-Smith, L., Kim, S., Calabretta, M.R., Haynes, R.F., 
& Kesteven, M.J. 2003, MNRAS, 339, 87 
Tenorio-Tagle, G., & Bodenheimer, P. 1988, ARA&A, 26, 145 
Thilker, D.A., Braun, R., & Walterbos, R.A.M. 1998, A&A, 
332, 429 

Walter, F., & Brinks, E. 1999, AJ, 118, 273 (WB99) 

Weaver, R., McCray, R., & Castor, J. 1977, ApJ, 218, 377 
Weaver, H., & Williams, D.R.W. 1973, A&AS, 8, 1 
Wouterloot, J.G.A., Brand, J., Burton, W.B., & Kwee, K.K. 
1990, A&A, 230, 21 


Appendix A: A comparison with other lists 

A.l. Heiles 

Heiles’ list of Hi shells in the Milky Way (Heiles 1979) 
contains fourteen shells in the 2nd Galactic quadrant iden¬ 
tified in the Hi survey of Weaver & Williams (1973). A 
visual inspection of the LDS shows that nine of the four¬ 
teen shells in the Heiles list are clearly and easily seen. 
Three shells are formed by partial walls and there are no 
well-defined holes. Two shells are completely open to the 
halo. 


Seven of the nine Heiles structures were re-identified 
by our method. The remaining two were found in veloc¬ 
ity channels but were later rejected because of unsuitable 
shapes of spectra. The derived properties of these nine 
structures (seven shells and two rejected candidates) are 
summarized in Table EH together with published values. 
The table shows a good correspondence both in angular 
dimensions and velocity extent. There are two cases where 
previous and newly found dimensions differ: GS091-04-69 
and GS128-I-01-105. GS091-04-69 is a shell which has in¬ 
ner walls — it might consist of several smaller shells — 
and we do not identify all of them as parts of the struc¬ 
ture. The other shell, GS128-I-01-105, is found to have a 
larger dimension in I and a larger range in radial velocity 
than was previously thought (see Fig. EH- left column 
and also Table EH- Judging by eye we think that in this 
case our derived dimensions are better. 

Three shells with partial walls were not identified as 
shells with our automatic method. We doubt the existence 
of one of them, GS123-I-07-127, as there is only a scarce 
group of fragments at the indicated position. The second 
shell, GS103-I-05-137, is very fragmented and we do not 
have any identification at this position, but we find a shell 
at these /, b coordinates in an adjacent velocity interval. 
According to Heiles, this shell is seen at velocities lower 
than —123 kms“^ and we see a structure in an interval of 
(-123,-100) kms“^. Very probably both identifications, 
Heiles’ and ours, are part of the same structure. The third 
structure from the list, GS130-I-00-I-15 is a large (35° in 
1) shell formed by a clearly observed wall, however this 
wall does not encircle the whole structure. The hole (the 
region of minimum in the temperature distribution) is not 
pronounced in this case. 

The two open structures (GS117-07-67 and GS129-05- 
91) are formed by incomplete arcs and therefore cannot 
be found by our algorithm. 


name 

Heiles 
deg^ kms”^ 

this 

deg^ 

paper 

kms-^ 

GS090-k02-115 

4x4 

X 

16 

3 

X 

4 

X 

12 

GS091-04-69 

9 X 10 

X 

28 

3 

X 

5 

X 

33 

GS091-k02-101 

4x3 

X 

12 

4 

X 

4 

X 

33 

GS108-04-23 

5 X 11 

X 

24 

7 X 14 

X 

20 

GS128-k01-105 

7 X 

6x4 

11 

X 

5 

X 

16 

GS148-01-H5 

4 X 

4x4 

3 

X 

5 

X 

10 

GS152-04-41 

4x4 

X 

12 

11 

X 

9 

X 

12 

GS095-k04-113 

GS139-03-69 

10 X 5 X 20 
18 X 10 X 28 

10 X 6 X 21 
13 X 6 X 22 


Table A.l. A comparison of previously known shells 
(Heiles 1979) and our identifications. The first column 
gives the names of the structures, the second one gives 
dimensions as published by Heiles, the third column gives 
our dimensions (AZ, Ab, Av). 
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Fig. A.l. Re-identification of the shell GS 128-1-01-105 discovered by Heiles (1979) - left column, and a newly discovered 
shell GS145-I-02-048 (Table 1) - right column. The upper row shows plots, the bottom row shows diagrams. 


A.2. Hu 

Hu (1981) made a list of high latitude shells detected in the 
Hi survey of Heiles and Habing (1974). The high latitude 
shells in Hu’s list are more difficult to see in LDS data; 
about half of them can barely be distinguished by eye or 
are simply not seen at all (this is connected to the method 
used for Hu’s identifications). The Hu list contains twelve 
shells in the second quadrant. Our automatic procedure is 
able to recover 50 % of them. Some shells from Hu’s list are 
too small (No. 26a) or fragmented (No. 26). Others (No. 
25 and 28) are not connected to any well-pronounced HI 
depressions. 

A.3. CGPS 

Our identification algorithm was tested on two fields of 
the CGPS, where two structures were found in the auto¬ 


matic search by Mashchenko & St.-Louis (2002). In the 
two data subsets, which were smoothed to a lower angu¬ 
lar resolution, we were able to recover structures Sh2-203 
and Sh2-187 with almost identical positions, dimensions 
and expansion velocities. 

A.4. Summary 

We identified with an automatic algorithm searching in 
the 2nd galactic quadrant in LDS 50% of the structures 
identified as Hi shells in a search in the Weaver and 
Williams (1974) survey performed by eye by Heiles (1979), 
and 50% of the high latitude shells discovered with a fil¬ 
tering method by Hu (1981) in the HI survey of Heiles 
and Habing (1974). We did not find shells formed only 
by incomplete arcs, or not connected to a minimum in 
lb maps. We found many more (nearly 300) shells than 
those contained in the Heiles and Hu’s lists. Many of these 
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newly identified shells are well defined (i.e. with a regu¬ 
lar shape and good contrast against the background) and 
have properties (dimensions, expansion velocities) com¬ 
parable to those previously found. An example of such a 
structure at the position: (Zc, be, wlsr) = (145.4, 2.1, -48.3) 
and {Al, Ab, Av) = (6.0, 4.5, 10.3) is shown in Fig. lA.ll 
(right column). Thus, it is surprising that they have not 
been previously discovered. 


